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ABSTRACT 

A variety of turbulence models, including five second-order closures and four two- 
equation models, are tested for the problem of homogeneous turbulent shear flow in a 
rotating frame. The model predictions for the time evolution of the turbulent kinetic 
energy and dissipation rate, as well as those for the equilibrium states, are compared with 
the results of physical and numerical experiments. Most of the two-equation models 
predict the same results for all rotation rates in which there is an exponential time 
growth of the turbulent kinetic energy and dissipation rate. The second-order closures 
are qualitatively superior since, consistent with physical and numerical experiments, 
they only predict this type of unstable flow for intermediate rotation rates in the range 
—0.1 < Cl/S < 0.6. For rotation rates outside this range, there is an exchange of stabilities 
with a solution whose kinetic energy and dissipation rate decay with time. Although the 
second-order closures are superior to the two-equation models, there are still problems 
with the quantitative accuracy of their predictions. 


* This research was supported by the National Aeronautics and Space Administration under NASA Contract No. NAS1-18605 
and was carried out while the author was in residence at ICASE, NASA Langley Research Center, Hampton, VA 23665. 
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1. Introduction 


The problem of homogeneous turbulent shear flow in a rotating frame constitutes a 
non-trivial test of turbulence models since it incorporates arbitrary combinations of shear 
and rotation which can have either a stabilizing or destabilizing effect. Since the flow 
is statistically homogeneous, two-equation turbulence models and second-order closure 
models give rise to an initial value problem for a coupled set of nonlinear ordinary 
differential equations. This allows for a fairly detailed mathematical analysis of the 
models based on dynamical systems theory and for the generation of inexpensive numerical 
solutions which have a high degree of accuracy. When these features are combined with 
the wealth of published results from physical and numerical experiments, a powerful 
tool for the testing and screening of turbulence models is obtained. This establishes the 
motivation for the present paper which constitutes an extension of the work of Speziale 
and Mac Giolla Mhuiris (1988). 

In this paper, a variety of two-equation models and second-order closure models will 
be tested for the problem of homogeneous turbulent shear flow in a rotating frame. The 
models to be considered are as follows: 

(i) the standard K — e model 

(ii) the nonlinear K — e model 

(iii) the RNG K — e model 

(iv) the Bardina K — e model 

(v) the Launder, Reece and Rodi model 

(vi) the Rotta-Kolmogorov model 

(vii) the Shih-Lumley model 

(viii) the Fu, Launder and Tselepidakis model 

(ix) the RNG second-order closure model 

This represents a broad collection of both two-equation models and second-order 
closures which includes some of the older models as well as some very recent ones. 
Comparisons will be made with the physical experiments of Tavoularis and Corrsin (1981) 
and the large-eddy simulations of Bardina, Ferziger and Reynolds (1983). It will be shown 
that several of the models yield reasonably good predictions for pure shear. However, the 
predictions of all of the models degrade considerably for appreciable values of the ratio of 
the rotation rate to the shear rate, f l/S. It will be shown that the two-equation models 
erroneously predict unstable flow (where the turbulent kinetic energy and dissipation rate 
grow exponentially in time) for all values of 0/5. The second-order closures will be 
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demonstrated to be qualitatively superior to the two-equation models since they predict 
unstable flow only for an intermediate range of fl/S; an exchange of stabilities occurs for 
other values of fl/S so that the turbulence decays consistent with the results of physical 
and numerical experiments. Although the second-order closures are able to account for the 
stabilizing or destabilizing effect of rotations on shear flow, the quantitative accuracy of 
their predictions is rather disappointing. The most surprising finding in this regard is that 
most of the newer models do not yield improved predictions over the older second-order 
closure models which have a considerably simpler structure. 

2. The Turbulence Models to be Tested 

Three of the two-equation models (i.e., the standard K — e model, the RNG K — e 
model, and the Bardina K — e model) are of the eddy viscosity type. For these models, 
the Reynolds stress tensor is of the form (see Hanjalic and Launder 1972) 


where 



is the mean rate of strain tensor, V{ is the mean velocity, K is the turbulent kinetic energy, 
and £ is the turbulent dissipation rate. The constant is typically taken to be 0.09 (in 
the RNG model, = 0.0837). In the nonlinear K — e model of Speziale (1987), a 
quadratic correction is made to equation (1). For flows in rotating frames, the Reynolds 
stress tensor for the nonlinear K — e model is given by (see Speziale 1989) 


n 3 = —KSij + 2 - 4 C D C^($ij 
+ Sik^kj “ mnSmn$ij + 2 Wikhj + 


K 3 /* 


where Cp = 1.68 and 


(3) 


o 



as, 


ij 


+ v. v Sij u ik S k j u>j k S k i 


u!{j + e m ji^lrn 




(4) 

(5) 

( 6 ) 


are the frame-indifferent Jaumann derivative, the intrinsic mean vorticity, and the local 
mean vorticity. In equation (5), flj is the rotation rate of the reference frame relative to 


2 


an inertial framing and e^-jt is the permutation tensor. 

The Reynolds stress models (l)-(3) must be solved in conjunction with the transport 
equation for the turbulent kinetic energy 

k -^w i - < 7 > 

(which is exact since the turbulence is homogeneous) and a modeled transport equation 
for the turbulent dissipation rate e. This modeled transport equation for e is of the general 
form (see Hanjalic and Launder 1972) 


■ r i - dVi 
e - c ^ K T a gx . 



For the standard K — s model and the nonlinear K — e model, 


( 8 ) 


C e i = 1-45, C e2 = 1-90 


( 9 ) 


whereas for the RNG K — e model of Yakhot and Orszag (1986) 

Cel = 1.063, C e 2 = 1.72 (10) 

In the Bardina K — e model (see Bardina 1988 and Bardina, Ferziger and Reynolds 1983), 
C e i an d C £ 2 are functions of W^j given by 

IS /l \ 1/2 

Cei = 1-50 - 0.015- UWijWiA (11) 

is / 1 \ 1/2 

Ce2 = 1.83 + 0.15- 1-WijWij J (12) 

The second-order closure models are based on the Reynolds stress transport equation, 
which for rotating homogeneous turbulent flows takes the form 

dvj dvi 

Tij ~~ Tik d^~ Tjk d^ k ~ 

where 11^ and £ij are the pressure-strain and dissipation rate correlations, respectively 
(again, ej j k is the permutation tensor). In the simplified form of the Launder, Reece and 
Rodi (1975) model that will be considered in this study, the models for II jj and £{j are 
as follows: 


{ T ik Q mkj^-m + T jk e 'rriki^Tn) II ij + £ ij (13) 
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( 14 ) 


£ ij = ( 15 ) 

where C\ and C 2 are dimensionless constants that are usually taken to be 1.8 and 0.6, 

respectively. Equations (13)-(15) are solved with the modeled dissipation rate equation 
(8) where C £ \ = 1.45 and C £ 2 = 1.90. 

In the Rotta-Kolmogorov model (see Mellor and Herring 1973) and the RNG second- 
order closure (Yakhot and Orszag 1988), the models for II jj and £{j are given by 

n y = Cl- | (ry + | Ml) + CiKSij (16) 

2 

£ ij (1^) 


where C\ = 3.21 and C 2 — 0.224 for the Rotta-Kolmogorov model and C\ = 1.59 and 
C 2 = 0.267 for the RNG second-order closure of Yakhot and Orszag (1988). Both the 
Rotta-Kolmogorov model and the RNG second-order closure are supplemented by the 
same modeled transport equation for e given in (8); in the Rotta-Kolmogorov model, 
C £ 1 = 1.8 and C £ 2 = 2.0 whereas in the RNG second-order closure, C £ \ = 1.063 and 
C £ 2 = 1.72. It should be noted that the Launder, Reece and Rodi model, the RNG 
model, and the Rotta-Kolmogorov model have one characteristic in common: the modeled 
pressure-strain correlation is linear in the anisotropy tensor 




(18) 


The Shih-Lumley model (1985) was constructed to satisfy realizability and incorporates 
more complex nonlinear anisotropic effects. The models for 11^- and take the form: 
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where 


+ 


/ 2 16 \ 
\3 + T as ) 


( dvk N 

r ik l ^ e mjk^rn 


+ r jk 



~^ e mik^rnj ~ ~^ T mnSmn^ij 


+ 5^{[ Tifc (a^ + 6mi ^ m ) 


+ Tjk 

(1^+e^iO™) 

T kl ^ T il r jmSlrn^ 

(19) 


£ . . _ -£$■■ 
£ ij ~ 3 eo »j 


(20) 

-ro( 

[l + 0.8 Vfj , 

F = 1 + 911 + 27III 

(21) 

i, 

~2 bij 

biji III = 

^jk^ki 

(22) 


F {80.1 In [1 + 62.4( — II + 2.3III)]} 

(23) 


The turbulent dissipation rate e is obtained from the modeled transport equation (8) with 
C s i and C £ 2 taken as 

Cgi = 1.44, C e2 = l + 0.49 [1 - 0.33 ln(l - 5511)] (24) 
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The Fu, Launder and Tselepidakis model (1987) is one of the newest models to be 
developed by Launder and his coworkers at UMIST. This model satisfies realizability and 
in the limit of small anisotropies reduces to the simplified version of the Launder, Reece 
and Rodi model discussed earlier. Its representation for II^- and £{j is as follows: 
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( 25 ) 


+ r{mi{T ik Wj k + Tj k W ik ) 

— 24b m ib n j {T mk W nk + T nk W mk )^\ 

£ ij = 6 ij - (l - Vf) Tij (26) 

where 

C{ = -6011, 7 = 0.6, r = 0.7 (27) 

The dissipation rate e is obtained from the standard modeled dissipation rate equation (8) 
with C e i = 1.45 and C £ 2 = 1.90. 

For homogeneous turbulent shear flow in a rotating frame (see figure 1), the mean 
velocity gradient tensor and the rotation rate of the reference frame are given by 


dvj 

dXj 


( 0 s 0\ 

0 0 0 , fii = (0, 0, f 2) 

\0 0 0 / 


(28) 


respectively. The nonlinear ordinary differential equations associated with each model 
are obtained in a straightforward fashion by substituting (28) into the tensor forms given 
above. These equations were solved using a Runge-Kutta-Fehlburg numerical integration 
scheme subject to the initial conditions 


K = K 0 , 


E £ 0 , T^j ^K 0 6ij 


(29) 


at time t = 0, which corresponds to an initially isotropic turbulence. Since the transport 
equations for r ij and e can be nondimensionalized into an equivalent set of transport 
equations for bij and e/SK, it follows that the solutions will only depend on the initial 
conditions through the dimensionless parameter e 0 /SK 0 (see Speziale and Mac Giolla 
Mhuiris 1988). 

3. Discussion of the Results 

The results obtained for the two-equation models will be discussed first. In figure 2, 
the time evolution of the turbulent kinetic energy K* = K/K 0 is shown as a function 
of dimensionless time t* = St for three rotation rates: Q/S = 0, £l/S = 0.25, and 
f l/S = 0.5. The predictions of the standard and nonlinear K — e model, the Bardina 
K — e model, and the RNG K — e model are compared with the results obtained from the 
large-eddy simulations of Bardina, Ferziger and Reynolds (1983) for an initial condition 
of e 0 /SK 0 = 0.296. It is rather disturbing to note that the newer models (i.e., the RNG 
and Bardina K — e models) yield substantially worse predictions than the standard K — e 
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model for pure shear as shown in figure 2a. While the RNG K — e model does well 
for the case where f l/S = 0.25 (see figure 2b), it appears that this is no more than a 
coincidence since the model erroneously predicts the same results for all rotation rates 
and considerably overpredicts the growth rate of the turbulent kinetic energy for the 
Cl/S = 0 and Cl/S = 0.5 cases (see figures 2a and 2c). The Bardina K — e model does 
yield results that depend on Cl/S\ however, the trends are in conflict with the large-eddy 
simulations. It predicts that the most energetic case occurs at Cl/S = 0 whereas large- 
eddy simulations indicate that the Cl/S = 0.25 case is the most energetic. Furthermore, 
in contrast to physical and numerical experiments, the Bardina K — e model does not 
predict a flow restabilization (where K and e — * 0 as t — ► oo) for all values of Cl/S that 
are sufficiently far outside of the range 0<fi/5<0.5. 

The equilibrium states of each of the two-equation models are compared in Table 1 with 
the physical experiments of Tavoularis and Corrsin (1981) and the large-eddy simulations 
of Bardina, Ferziger and Reynolds (1983). It should be noted that since the LES results 
were not defiltered, they tend to considerably overpredict the anisotropy tensor and should 
only be considered for the trends that they exhibit. All the equilibrium states were 
computed with C ^ « 0.09. Even though improved results are obtained when is 
considered to be a function of the ratio of production to dissipation as suggested by Rodi 
(1972) (which yields (C^)oo = 0.055 as used by Speziale and Mac Giolla Mhuiris 1988), 
in this study we will only show the equilibrium results corresponding to C ^ = 0.09, since 
this value is almost always adopted in applications. Several observations concerning these 
results are noteworthy: (a) both the standard and nonlinear K — s models yield the same 
values for {SK/s)oq and (612)00 which are barely in the range of the physical experiments 
for pure shear, (b) only the nonlinear K — e model is able to reproduce the trends of 
the large-eddy simulations for the normal components (6n)oo and (622)00) and (c) the 
RNG model drastically overpredicts the values of (612)00 and (SK/e) qq in pure shear 
which accounts for the overly large growth rate of the turbulent kinetic energy shown in 
figure 2a. 

The equilibrium diagrams of (e/SK) 0 o are shown as a function of Cl/S in figure 3 for 
the standard, nonlinear, RNG and Bardina K - e models. It is clear from this figure 
that all of the two-equation models considered in this study are incapable of properly 
accounting for the effects of rotations on turbulent shear flow. The standard, nonlinear, 
and RNG K — e models erroneously predict the same results for K and e independent of 
the rotation rate. While the Bardina K — e model yields equilibrium states that depend 
on Cl/S, they appear to be unphysical. For example, there are no stable equilibrium 
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solutions for rotation rates outside of the range —0.9 < Cl/S < 1.9; thus, the model does 
not predict the flow restabilization that is known to occur for most Cl/S outside of the 
range 0 < Cl/S < 0.5. 

The time evolution of the turbulent kinetic energy predicted by the various second- 
order closure models is shown in figure 4 for the same three cases discussed above. Three 
important conclusions can be drawn from these results: (a) the Shih-Lumley model and 
Rotta-Kolmogorov model predict too small a growth rate for pure shear whereas that 
predicted by the RNG model is too large in the initial stages (see figure 4a), (b) all of 
the models yield rather bad predictions for the most energetic case of Cl/S = 0.25 (see 
figure 4b), and (c) the Launder, Reece and Rodi model, the Shih-Lumley model and the 
Fu, Launder and Tselepidakis model erroneously predict stable flow at Cl/S = 0.5 whereas 
the RNG second-order closure has much too strong a growth rate (see figure 4c). 

The equilibrium values predicted by each second-order closure is compared in Table 2 
with the physical experiments of Tavoularis and Corrsin (1981) and the large-eddy 
simulations of Bardina, Ferziger and Reynolds (1983). The performance of these models 
can be summarized as follows: (a) the Launder, Reece and Rodi model does reasonably 
well in predicting the equilibrium values for pure shear but is defective in its prediction of 
a premature restabilization at Cl/S « 0.39, (b) the Rotta-Kolmogorov model erroneously 
predicts the opposite extreme wherein the Cl/S = 0 and Cl/S = 0.5 case are equally 
energetic since they have the same values for (612)00 and (SK/e) 00, and (c) the Fu, 
Launder and Tselepidakis model does the best in its prediction of the equilibrium states 
for pure shear; however, it also erroneously predicts a premature flow restabilization 
at Cl/S « 0.39. It should be noted that the erroneous prediction of (612)00 = 0 and 
(SK/e)<x) = 00 made by the RNG second-order closure is directly tied to its rather low 
value of C £ \ = 1.063. This model could become more competitive with the other second- 
order closures if C e \ was of the order of 30% larger. 

A characteristic equilibrium diagram for the second-order closures is shown in figure 5 
(note that all the models possess the same topological structure for (e/SK) 00 taken as a 
function of Cl/S). The models have two equilibrium solutions: one where (e/SK) <30 7^ 0 
(the elliptic region AB) for which e and K grow exponentially with time at the same rate, 
and one where (e/SK)oo = 0 for which e and K predominantly decay with time. Hence, 
the second-order closures properly predict unstable flow only for intermediate rotation 
rates which, depending on the model, can lie anywhere in the ranges —0.11 < Cl/S < 0.39 
to —0.09 < Cl/S < 0.61. In this intermediate band of Cl/S, the two solutions undergo 
an exchange of stabilities which qualitatively mimics the shear instability. Furthermore, 
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the second-order closures predict a universal equilibrium (i.e., one that attracts all initial 
conditions) in the unstable regime consistent with the experiments of Tavoularis and 
Corrsin (1981). 

4. Concluding Remarks 

It has been demonstrated in this study that the existing two-equation turbulence 
models of the K — s type do not properly account for the stabilizing or destabilizing 
effect of rotations on turbulent shear flow. Furthermore, the quantitative predictions of 
many of these two-equation models are rather poor for the simple case of pure shear flow. 
It appears that the nonlinear K — e model yields the best results among the two-equation 
models, however it is still highly deficient. The second-order closures are far superior 
in a qualitative sense since they do properly predict that there is unstable flow only 
for intermediate rotation rates that are not too far removed from the unstable regime 
of 0 < f l/S < 0.5 predicted by linear stability theory. However, two major problems 
remain with the quantitative accuracy of the results of the second-order closures: (a) the 
prediction of a premature flow restabilization at 0,/S « 0.39, or (b) the prediction of poor 
equilibrium values for the case of pure shear, f 't/s = 0. Since it is important to predict 
simple shear flows accurately in engineering applications, it is our opinion that (of the 
models considered in this study) the Fu, Launder and Tselepidakis (1987) model is on 
balance the best. However, like the other existing second-order closures, it does not yield 
accurate results when there is both shear and rotation. The development of significantly 
improved models based on invariance arguments and a dynamical systems approach will 
be the subject of a future paper. 
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Equilibrium 

Values 

Standard 
K-er Model 

Nonlinear 
K-c Model 

RNG 

K-e Model 

Bardina 
K-e Model 

Experiments 

Large-Eddy 

Simulations 

o 

o 

II 

cT 

b n 

0 

0.353 

0 


0.201 

0.30 

b 2 2 

0 

-0.252 

0 


-0.147 


bn 

-0.212 

-0.212 

-0.489 

HEE2HI 

-0.150 


WESBtk 

4.71 

4.71 

11.70 

5.49 

6.08 

- 

tl/S = 0.25 

bn 

0 


0 

0 

- 

0.06 

b 2 2 

0 


0 


- 

0.05 

bn 

-0.212 


-6.489 

■ 

- 

-0.35 


4.71 

4.71 

11.70 

4.84 

- 

- 

n/s = 0.5 


0 

-0.252 

0 

0 

- 

-0.26 

6 2 2 

0 

0.353 

mm 

0 

- 

0.25 

bu 

-0.212 

-0.212 

-0.489 

-0.194 

- 

-0.10 

SKje ) 

4.71 

4.71 

11.70 

4.30 

- 

- 


Table 1. Equilibrium results for homogeneous turbulent shear flow in a rotating frame. Comparison of the predic- 
tions of the two-equation models with the Large- Eddy Simulations of Bardina, Ferziger, and Reynolds 
(1983) and the experiments of Tavoularis and Corrsin (1981). 



Equilibrium 

Values 

LRR 

Model 

RK 

Model 

SL 

Model 

■ 

RNG 

Model 

Experiments 

Large-Eddy 

Simulations 

^3 

'co 

II 

0 

6n 


0.242 

0.114 

0.214 

0.533 

0.201 

0.30 

&22 


-0.121 

-0.116 


-0.267 

-0.147 

-0.26 

612 


-0.169 

-0.122 


0 

-0.150 

-0.15 



3.71 

6.93 

7.28 

00 

6.08 

- 

n/s = 0.25 

bn 

-0.059 


-0.070 

-0.134 

0.159 

- 

0.06 

b 22 

0.155 


0.067 

0.200 

HI 

- 

0.05 

bn 

-0.208 


EsEESI 

-0.160 

-0.449 

- 

-0.35 


4.83 


IUJJ1 

6.54 

12.70 

- 

- 

10 
d 

11 

cf 

*.1 

-0.120 

-0.121 

-0.123 

-0.1511 

-0.267 

- 

-0.26 

622 

0.160 

0.242 

0.120 

0.061 

0.533 

- 

0.25 

612 

0 

-0.169 

0 

0 

0 

- 

-0.10 

SKje 

00 

3.71 

00 

00 

00 

- 

- 


Table 2. Equilibrium results for homogeneous turbulent shear flow in a rotating frame. Comparison of the pre- 
dictions of the second-order closure models with the Large-Eddy Simulations of Baxdina, Ferziger, and 
Reynolds (1983) and the experiment of Tavoularis and Corrsin (1981). (LRR = Launder, Reece and 
Rodi, RK = Rottar Kolmogorov, SL — Shih-Lumley, FLT = Fu, Launder and Tselepidakis). 
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ison of the predictions of the two-equation models with the Large-Eddy Simulations of 
Bardina, et al. (1983): (a) 0/5 = 0, (b) 0/5 = 0.25, and (c) 0/5 = 0.5. 

Figure 3. Equilibrium diagram for the two-equation turbulence models in rotating shear 
flow. 

Figure 4. Time evolution of the turbulent kinetic energy in rotating shear flow. Compari- 
son of the predictions of the second-order closure models with the Large-Eddy Simulations 
of Bardina, et al. (1983): (a) 0/5 = 0, (b) 0/5 = 0.25, and (c) 0/5 = 0.5. 

Figure 5. Equilibrium diagram for the second-order closure models in rotating shear flow. 
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